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Abstract 

The mass-constraining variable M 2 , a (1 + 3)-dimensional natural successor of extremely popular 
Mt 2 i possesses an array of rich features having the ability to use on-shell mass constraints in 
semi-invisible production at a hadron collider. In this work, we investigate the consequence of 
applying a heavy resonance mass-shell constraint in the context of a semi-invisible antler decay 
topology produced at the LHC. Our proposed variable, under additional constraint, develops a 
new kink solution at the true masses. This enables one to determine the invisible particle mass 
simultaneously with the parent particle mass from these events. We analyze in a way to measure 
this kink optimally, exploring the origin and the properties of such interesting characteristics. We 
also study the event reconstruction capability inferred from this new variable and find that the 
resulting momenta are unique and well correlated with true invisible particle momenta. 
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I. INTRODUCTION 


The Large Hadron Collider (LHC), after the successful completion of its first run and 
the remarkable achievement of the Higgs boson discovery BE], has already entered into 
its second phase. Upgraded with higher energy and luminosity, the main physics goal is 
to explore the multi-TeV scale associated with the physics beyond the Standard Model 
(BSM). Although the LHC has not reported any clinching evidence for new physics so far, 
expectations are running high for possible new physics signals in soon unless such signatures 
are already hidden inside the LHC data. Any scenario with a positive outcome essentially 
demands the measurements of the mass, coupling and spin of new BSM particles. However, 
this is going to be complex since many of the very likely scenarios with a wide class of BSM 
models have embraced the concept of thermal relic dark matter (DM) as some stable exotic 
member within them. Since these massive DM particles are colorless, electrically neutral 
and weakly interacting, once produced in the collider, they do not leave any trace at the 
detector. Hence, one needs to rely on the experimentally challenging signature of missing 
transverse momenta /V from the imbalance of total transverse momentum, accounting for 
all visible decay products in each event from an already adverse jetty environment of the 
hadron collider. Moreover, the DM in many models is expected to produce in pairs because 
of its stabilizing symmetry,^ commonly associated with the Z 2 parity, which makes the event 
reconstruction even more challenging. 

Keeping an eye on the above encumbrances, many ideas have already been developed for 
the determination of mass and spin. For some recent reviews, see Refs. for different mass 
measurement techniques. To give a brief and sketchy classification for the mass measurement 
techniques, parts of which are also crucial for our study, one looks into the following: (a) 
In the end-point method [TOl - fiG] . the end point of the independently possible invariant 
mass constructed out of all combinations of visible particles is related to the unknown 
masses involved in a decay chain. Thus, a sufficient number of end-point measurements 
are needed to pin down all the unknown masses, which is only possible for a longer decay 
chain, (b) The polynomial method [I7H2T] tries to solve for all unknowns (unknown masses 

^ There can be a DM stabilizing symmetry other than Z 2 which allows more than one DM particle per 
vertex and kinematic variables [sm were used to distinguish between the DM stabilizing symmetries, but 
we stick to the popular Z 2 parity for simplicity and extendability. 
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and invisible particle fonr momenta) in the event, using all available constraints. Again 
all unknown masses can be determined only if one has at least a three-step decay chain, 
(c) The transverse mass methods are based on transverse mass variables dehned by the 
minimum parent mass consistent with minimal kinematic constraints in the event, for a 
given trial mass of the invisible particle. Here constraints consist of the equality of parent 
mass, mass-shell relations and the missing transverse momenta from undetected invisible 
particles. Several interesting features of these variables received lots of attention recently. 
The transverse mass variables were further extended with many variants, such as, Mt 2 & 
150] , subsystem Mt 2 [31] , Mt 21 . and its sister Mt 2 \\ [32] , Mct 2 [331133] , asymmetric Mt 2 [331 
[36] . contransverse mass Mct [37113^ . Mct± and its sister Mct\\ [3D]- In Ref. [H] I32] it 
is shown that the (1 -|- 3)-dimensional generalization of Mt 2 can be useful, where one can 
apply the equality of parent mass and/or equality of relative^ mass constraints to improve the 
number of events at the end point, (d) Inclusive variables are the ones whose construction 
do not take into account the production mechanism or topology of the new particles and 
only depend on the hnal state visible particles momenta and missing transverse momenta 
in the event. There are many inclusive variables, such as Ht [33], total visible invariant 
mass M [33], -^T, effective mass Mg// [10], and its sisters [3SH3E]- They are useful for 
extracting the mass scale of new physics in a model independent way. While each method 
is more relevant for a particular outlook, the transverse mass variable and its extensions 
turned out to have the potential for precise measurement of the masses in short single-step 
decay chains involving invisible particles. Consequently, for our present analysis, we further 
develop and discuss this method in the next section. 

The present work is inspired from another class of prescriptions known as the cusp 
method |39ll52] . which can determine both the parent and invisible particle (or DM) masses 
simultaneously from single-step decay chains. However, the primary assumption is that both 
the parents are produced from an on-shell heavy resonance, whose mass is known already. 
This method uses the measured momenta of visible particles and constructs possible vari¬ 
ables, among which the invariant mass and magnitude of the individual transverse momenta 
of the visible particle can exhibit non-smooth structure or cusp and end point in their distri¬ 
butions. The location of this cusp and end point are function of the unknown masses, so the 

^ Relative is defined by any resonance other than parent and daughter realised in some particular topology. 
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parent and danghter mass can be determined simnltaneously from these two measurements. 
It is argued that the longitudinal boost invariant variables, such as the individual transverse 
momenta and the invariant mass of visible particles, are minimally affected by the unknown 
transverse boost, which are sufficient to pin down both the unknown masses. This analysis 
was extended further in the context of the linear collider |52]. 

In this paper, we propose a complementary procedure to measure both the parent mass 
and daughter mass simultaneously using the M 2 variable. We also assume that the heavy 
resonance mass is known, and this information is embedded in the minimization of M 2 
resulting in a new constrained variable dubbed M 2 Cons- Further discussion on this variable 
is in Sec. |TT] with formal definitions and chronological motivations. Emboldened by the 
interesting characteristic of this variable which develops a kink in the distribution end-point 
maxima corresponding to the correct value of the unknown invisible mass, we propose a 
simultaneous measurement of both the masses by identifying the kink position. We also put 
forward a desirable way of recasting the kink by utilizing all available data, giving a robust 
outcome irrespective of other realistic effects. This M 2 Cons variable proved to be useful in 
all mass range, including the region known as the ‘large mass gap’ in the kinematic cusp 
method where the cusp may not be very sharp. A unique event reconstruction can also be 
realized using M 2 Cons- 

As a (l-|-3)-dimensional variable, M 2 has the capability to use all components of three 
momenta. Thus, the utilization of additional mass constraints is possible if such information 
is available. This brings additional advantages, contrary to the transverse mass variables, 
which by construction are not similarly capable. Here we argue that the newly developed 
M 2 Cons variable equips one to consider additional constraint from on-shell mass resonance, 
which in turn produces suitable kinematic constraints to get a new kink solution in antler 
topology. Importantly, this contribution towards the kink is not merely symbolic but is 
substantive, enabling better mass measurements. 

This paper is organized as follows. In Sec. [IT], we motivate the antler production topology 
at the hadron collider describing all available constraints associated with such topology. In 
this context, we also open up discussion on the Mt 2 transverse mass variable and its (1-1-3)- 
dimensional generalized sisters in the M 2 family. Following that, we introduce our new 
variable M 2 Cons- We further discuss the effect of an additional heavy resonance constraint. 


describing the basic features and benefits of this variable. Sec. Ill is assigned to describing 
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the kink solution coming from the M 2 Cons distribution end point as a function of the trial 
invisible mass which is yet unknown. Further comparison is made with the corresponding 
maximum achievable in the M 2 variable. Realizing that the number of events at the end 


point is very limited for M 2 Cons, we formulated an efficient way in Sec. IV to reconstruct 
the kink by recasting all available data. In Sec. |V| we show that M 2 Cons can be used for 
a unique event reconstruction, and the reconstructed momenta are well correlated with the 


true momenta of the invisible particle. Finally, in Sec. |V^ we summarize our main results 
and conclude. 


II. ANTLER TOPOLOGY AND CONSTRAINED VARIABLE WITH M2 

Antler topology, realized in different Standard Model (SM) processes and a variety of 
new physics models, is very common and widely explored. The SM Higgs boson decaying 
semi-invisibly through the W-boson,^ h —)■ IF + IF* -\- iu, or via r lepton, h —)■ 

r + r —)■ W*i'r + W*Ur, are some of the signihcant channels in the context of the SM Higgs 
search at the LHC. Similarly, in several BSM theories, the search strategy relies on antler 
production topology. Some of these include the heavy Higgs of supersymmetry (SUSY) 
decaying to the Z-boson and lightest supersymmetric particle (LSP) via neutralinos, H —)■ 
X 2 + X 2 |S3] and the SUSY extended Z' decaying to the lepton and LSP via 

the slepton, Z' —>■ £+ + £“ —>■ + [SUES]. Similarly, in a universal extra-dimensional 

model, second excitation states can decay to hrst excitation states, Z^‘^'> —)■ Pd) _|_ /^(i) 
£-^(1) £+^(1) j56[ |5^. The semi-invisible decay of doubly charged exotic scalars in many 

BSM scenarios can produce SM particles via W pairs, —)■ IF''' + IF’'' —)■ [58] . 

Moreover, the heavy Higgs or heavy Z' can also decay semi-invisibly to SM particles via ti 
pairs, H/Z' + bW^ + bW~ —)■ bi'^ug + In addition, antler topology can also 

be realized at the linear collider as hxed c.m. energy is equivalent to the heavy resonance 
produced at its rest frame before pair production and subsequent decay (for example, see 

H). 

Before starting our analysis, let us describe the basic setup and the notation. The rep- 

^ Note that except from the fact that h ^ W + W* probably signify most familiar SM antler channel, this 
off-shell production of IF is not pursued further in present analysis. 
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FIG. 1: Archetype of antler topology where G, a heavy resonance particle with mass mo produced 
at hadron collider, decays to two daughter particles Pi and P 2 through two-body decay, each 
of which subsequently decays to produce SM visible particle (Vi) and an invisible or dark matter 
particle {xi)- Momenta of these visible and invisible particles are assigned as pi and qi, respectively. 
mp and ruy- are masses of parent Pi and invisible particle Xi- This topology can be considered for 
SM Higgs production with subsequent decays into SM visible and massless neutrinos as invisible 
particles. On the other hand, in a BSM scenario, G (a parity-even state) can decay to produce 
(parity-odd states) Pi and Xi- 

resentative diagram for antler topology is shown in Fig. [T| where a parity-even"^ heavy 
resonance particle G (grandparent) with mass me decays to two parity-odd particles Pi and 
P 2 , each of which snbsequently decays to the Standard Model particle (V)) and an invisible 
or dark matter particle {Xi)- We assign the momenta to visible and DM particles in each 
side of the decay chain as Pi and qi, respectively. Moreover, we denote the masses of parents 
{Pi) and invisible danghters (xi) as mp and m^, respectively. The primary motivation of 
this analysis is to determine these nnknown parameters. Thongh we have shown a generic 
antler topology in Fig. in this analysis we are interested in the symmetric antler process 
motivated by the above examples. The symmetric antler inclndes same parent (Pi = P 2 ) 
and same danghter {xi = X 2 ) particles, or at least their masses are same, = mp^ = mp 
and = m^^ = m^. 

We wonld like to address this topology nsing the mass-constraining variable for the snb- 
^ Parity is pertinent only for the BSM processes having stable invisible exotic particles in the hnal state. 
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system represented by the gray shaded region shown in Fig. [Tj Let us start with the existing 
and popular transverse mass variable Mt2 before moving into the generalization and hnally 
extending to our new variable. Mt2 is dehned to have the potential to measure the masses 
of the BSM particles both in short or long decay chains, although its dominance and signih- 
cance is mostly grounded in its capability to handle the former case. The classic dehnition® 
of Mt2 is given by the larger value between two transverse masses ^ constructed from 
both sides of the decay chain and minimized over unknown invisible momenta satisfying the 
Pt constraints of that event. Mathematically, 


Mt2 = min 


QiT 

QiT=i^T} 


max{M!p\piT, qiT, 

2 = 1,2 

with the usual dehnition of the transverse mass for each decay chain. 


( 1 ) 


{mPY — + 2{Ej, 


\visii) T-imttfi) 

tjrj 


PiT-qiT) 


e: 


,vis(i) 




inv{i) 


( 2 ) 

( 3 ) 


{Ep^^^\piT) and qip are (l+ 2 )-dimensional transverse energy-momenta correspond¬ 

ing to the visible and the invisible decay products in the decay chain, respectively. Note 
that, in the dehnition of Mt2-, the minimization is done over all possible partitions of /V 
and the maximization of within the bracket ensures a closer shot towards the parent 
mass mp. By this dehnition, Mt2, calculated for each event, must be smaller than or equal 
to mp. 

The maximum quantity for any of these mass variables M,.. (such as, Mp2, M2 or M2Cons) 
over the available data set is 


\All events) 

Now should provide a very close estimate of mp. Moreover, in the scenario where 

the invisible particle mass is a priori unknown, e.g. dark matter models, {fhP) would 

still oher a useful correlation with the trial invisible mass a daughter mass hypothesis 

® Here, “T” and “2” in Mp 2 stands for the transverse projection and two parent particles, respectively, in 
the topology under examination. Reference [5] generalized and unified the concept of mass variables and 
set a preferred nomenclature according to the order of operations to rewrite the same variable as M 2 T 
within a general M 2 family. Notably, [S] also demonstrated the fact that the transverse projection can be 
done using not one, but three completely different schemes. 
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used as an input for the calculation of Mt2- One can possess only this partial information 
on the unknown parent and daughter masses unless, under some special circumstances, this 
correlation curve generates a kink feature exactly at the correct mass point. We will have a 
slightly elaborate discussion about the kink feature in Sec. |III[ 

Now to motivate the (1 + 3 )-dimensional generalization of previous dehnitions as in Eq. 
one readily notes that the Mt2 is not utilizing longitudinal components of the momenta and, 
thus, the available mass-shell constraints for a given topology. M2 is thus constructed | 9 ] 
out of the (1 -|- 3 )-dimensional momenta by removing all “T” in the dehnition of Eqs. [T]j^ 
(except that of the total missing transverse momentum constraint under the curly bracket, 
since the longitudinal part is not available in the context of the hadron collider). Now one 
can apply the on-shell mass constraints in the minimization of M2, and, depending on the 
constraints applied, different constrained classes of the M2 variable {e.g. M2xx, M2cx and 
^2cc) can be constructed; details about these variables can be found in Ref. [ 12 ]. Using 
similar notation, one can readily come up with the hrst two types of variables available from 
the subsystem considered in Fig. Here, M2cx is the (1 + 3 )-dimensional generalization of 
Mt2 with the equality of the parent mass constraint applied in the minimization. 


M2CX = 


mm 

qim 

1 <flT+<?2T=^T 

l(pi+gi)^=(p2+<?2)^ 


max{MW(pi, 

i=l,2 


with the (1 -|- 3 )-dimensional mass from each decay chain as 


( 5 ) 


- p, 4 ). 


( 6 ) 


The corresponding M2xx variable is simply perceived once the last constraint inside bracket 
is absent, just like the transverse mass case in Eq. It is straightforward to show [ 12 ] 


Mt2 — M 2 XX = M 2 
= M 2 cx- 


( 7 ) 

( 8 ) 


Also note that, in our example, there is one visible particle per decay chain in the hnal state. 
Hence, the Mt2 and other variables always come from a balanced conhguration irrespective 
of the choice of trial invisible mass. So once again the maximum Mtpf^ (or the maxima of 
other variables as in Eq. and Eq. can only give a constraint between parent and invisible 
particle mass. 




Now, by following the steps before the subsystem as in Fig. [T| one realizes that the parents 
(Pi, P2) are actually originated from a heavy resonance {G). In a BSM scenario, even-parity 
G can directly decay to SM observable particles and hence the mass, me, can in principle 
be measured. Here, before we move further, we assume that in our topology only this heavy 
resonance mass me is known. We are now in a position to develop a variable using this 
mass constraint, so that® 


^2Cons 


min 


Qit+Q2t=I/t 

ipi+P2+qi+q2)^=m^ 


max{MW(p.,g.,7u^.^ 

i=l,2 




(9) 


where (1 -|- 3 )-dimensional invariant mass is as in Eq. Additionally, the dependence 
on the unknown trial invisible mass m^ is shown explicitly. With this additional constraint, 
one expects a squeezed phase space affecting this new variable from that of M2; furthermore, 
we will soon realize that this effect is a little more far-reaching. Before we gradually move 
to demonstrate that, let us open the discussion with the consequences of this new variable 
in the invisible momenta space. 

The additional heavy resonance mass-shell constraint in the minimization (last condition 
inside bracket of Eq. constrains the invisible particle momenta, such that, the invariant 
mass of the parents is conhned to a narrow resonance of mass tuq. This phenomenon is 
true for each event. In Fig. the effect of this constraint is demonstrated for one event 
with an example where the trial invisible particle mass is considered smaller than the 
yet unknown true mass The region represented by the light temperature map color 
gradient is the maximum between two transverse masses M^\ as in Eq. [^before executing 
the minimization. This is shown with respect to the invisible momenta components, qi^ and 
qiy, by taking care of the missing transverse momenta constraints. Now the minimum of this 
quantity, which is nothing but the Mt2, is the minimum point in the color map displayed 
by the hlled circle •, and different contour lines are shown by dashed curves. Moving to our 
(l-|- 3 )-dimensional new variable with the heavy resonance mass-shell constraint in Eq.|^ once 
again after doing the similar exercise we get the solid contour curves superimposed in the 
same plot. Of course, we are no longer showing the color gradient as done in the transverse 


® One can also consider an additional constraint using the equality of parent mass (pi + qi)'^ = (p 2 + 92 )^ 
in M 2 Cons ■ Although this would further constrain the allowed invisible momentum space, it would finally 
choose the same minima. Hence, our arguments with this present example and analysis remain the same. 
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mass case. Transformation of dashed contonr lines into corresponding solid ones (same color 
represents the same valne of that contonrs) within the same region qnalitatively indicate the 
effect of this additional mass constraint. Minimnm of these solid contonrs represents the 
M 2 Cons) displayed by circle pins © in the same hgnre. Note that the longitndinal momenta 
components for the invisible pairs are eliminated in this demonstration by minimizing with 
the G mass-shell constraint. 

As noted, mo constraint restricts the invisible momenta making the region shrink as 
depicted by the dashed and solid contonr {e.g. following bine lines correspond to 100 GeV), 
the dashed line contonr does not satisfy the additional G mass-shell constraint while solid 
contonr does. The same is trne for all other lines also. The represented valnes of Mt 2 and 
M 2 Cons considered in this example are 75.1 and 98.5 GeV, respectively, for trial mass rh^ at 
10 GeV which is smaller than the trne invisible mass 100 GeV. The corresponding trne mass 
of parents and heavy resonance are 200.0 and 1000.0 GeV, respectively. The white dashed 
(solid) line represents the eqnality of the transverse-mass (mass) of parents and this eqnality 
line has moved towards higher valne becanse of the constraint. As a resnlt, one natnrally 
expects M 2 Cons > M 2 event by event. The red star ★ is the position of the trne transverse 
momenta of the invisible particle. Glearly, constraint brings the minima, M 2 Cons, closer to 
the trne momenta, and this can improve any effort to reconstrnct the invisible momenta. 
This featnre will be farther considered and discnssed in Sec. lYl 

One more remarkable featnre emerges at this point which will be capitalized in the next 
section. We already noted the event wise npward shift of valnes nnder the constraint. The 
next natnral qnestion in this context concerns the maximnm valne achievable by this mass 
variable and how it is related to the trial missing particle mass? The experimentally mea- 
snred maxima can deviate (downward) from the theoretical maximnm of the mass variable 
depending npon the accessible nnmber of events, and more importantly, the abnndance of 
events towards the end point of the distribntion. We postpone this issne for the time being 


and consider it again in Sec. |IV[ Now coming back to onr variable, it shonld not be snrprising 
that at the trne valne of the invisible particle mass (i.e. when = m^), the maximnm 
valne of the constraint variable coincides with that of the variable withont this con¬ 
straint, This is becanse both of these variables are derived for the same topology. 

On the other hand, at all other trial mass valnes, not only is the individnal (event-by-event) 
constraint qnantity larger, bnt also the maximnm of that constraint mass is the larger valne. 
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FIG. 2: The effect of the heavy resonance mass-shell constraint is demonstrated using the mass 
variables Mt 2 and M2Cons considering one antler event in the invisible momentum component 
space. The region represented by the color gradient is the maximum among two transverse masses 
coming from two decay chains. Corresponding contours are shown with dashed lines. Following 
Eq. the minimization of this quantity, the Mt 2 ■, is represented by the hlled circle #. In the same 
plot, the solid lines (of same colors) are delineating the corresponding contours for the (1 + 3)- 
dimensional new variable with heavy resonance mass-shell constraint as in Eq. Note that only 
the contour lines are shown in this case, not the color gradient as in the transverse mass case. 
The minimum of these solid contours is represented by the M 2 Cons, displayed by circle plus © in 
the same hgure. The G mass-shell constraint restricts the invisible momenta, making the region 
shrink, as depicted by the dashed and solid lines. The white dashed (solid) line represents the 
equality of transverse-mass (mass) of parents and this equality line is also moving towards higher 
values because of the constraint. The red star ★ is the position of true transverse momenta of 
invisible particle. The mass spectrum we choose is {me, mp, = (1000.0, 200.0, 100.0) in GeV 
and the trial invisible particle mass we took for this plot is 10.0 GeV 


To write in a compact form, 

M max 

f 


= M^°'^{m^) = mp, if 


2Cons['^x) 


( 10 ) 


X ~ ""x 

> M^^^{m^), if 

While this point is further discussed in the following section as a means of measuring the 
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FIG. 3: Normalized distributions of the mass variables are delineated using a toy model of antler 
topology with parents mass at 200 GeV. Both the M 2 (dark) and M 2 cons (green) distributions, 
considering the invisible particle mass at its true value (100 GeV), produce the end point at the 
correct parent mass. However, the heavy resonance constraint gives M 2 cons a higher value, resulting 
a larger number of events at the end point. 

unknown masses, here we illustrate it with one example distribution for the aforementioned 
M 2 variables, considering a toy process with the antler topology as shown in Fig. [T} For 
demonstration purposes, we choose a mass spectrum with {me, mp, = {1000, 200,100} 
in GeV, which is a relatively difficult region for the kinematic cusp method l50] known as 
the “large mass gap” region, where the cusp may not be very sharp, leading to large errors 
in the mass determination. The M 2 Cons variable can be effective for mass determination 
both in the large mass gap region as well as in other regions of phase space. In Fig. 
we have compared the normalized distributions for M 2 and M 2 Cons at the true mass of 
the invisible particle. Both the constrained (green histogram) and unconstrained (dark 
histogram) distributions share the same end point precisely at the parent mass, as argued 
earlier. However, from the distribution, one should also note the movement of the events 
towards the higher value under the heavy resonance constraint and, thus, expect a larger 
number of events at the end point. 
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III. MASS MEASUREMENTS WITH KINK 


Kink in mass measurement techniques is a widely acclaimed feature, first shown in the 
context of the Mt 2 variable. Let us continue from the brief discussion below Eq. where 
the distribution maximum (rh^) is observed to offer a useful relation correlating the 

parent mass with the trial value of unknown invisible mass It is shown [261 EH that 
simultaneously both parent mass mp and daughter mass can be determined by iden¬ 
tifying a kink (continuous but not differentiable) in this correlation curve, where the true 
mass point resides. It is worthy of attention that the {fh^) has two different functional 
forms before and after this kink, and they share the same value at the true mass point. This 
behavior stems from the fact that the visible system invariant mass of any (or both) decay 
chain(s) have to have a range of values; hence, there should be at least two visible particles 
per decay chain. Consequently, experimentally simpler single-step decay chain topology is 
deprived of such advantage. The above feature was shown where the system does not have 
any recoil from initial state radiation (ISR) or upstream transverse momenta (UTM). But 
the presence of ISR is inevitable during the production at any hadron collider. It is sub¬ 
sequently revealed [2H1E91I3I] that kink can also arise from topology having a single-step 
decay chain on both sides, but there should be recoil to the system which may come from 
ISR or UTM. Both scenarios can naturally arise in subsystem context from a longer decay 
chain. However, sizable kink resolution only comes from the events with very high recoil 
Pp, essentially with very low statistics. 

In the last section, we dehne the constrained variable M 2 Cons using the heavy resonance 
on-shell constraint in the minimization of M 2 . Analogous to the previous discussion,^ the 
maximum of this variable also exhibits different dependence at either side from true 


mass, as a function of trial invisible particle mass m^. Following the Eq. 10, one can obtain 
the kink structure exactly at the true mass. However, the source for the appearance of 
this kink is attributed to the heavy resonance mass-shell constraint in the minimization.® 


^ At this point, we would like to make it clear that the effect of ISR/UTM is not considered in this present 
analysis. This study shows a new kink solution due to the kinematic constraint coming from on-shell 
mass resonance in antler events. If one consider such events associated with ISR, that may marginally 
contribute strengthening the case over already strong kink solution as demonstrated. 

® One can argue that the heavy resonance constraint works in the same spirit of ‘relative’ constraint as 
defined in M 2 class of variables [42]. In fact, in a non-antler scenario the M 2 variables under usual relative 
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FIG. 4: Upper plot depicts the behavior of the upper end points for both the constrained 
unconstrained variables with respect to the trial invisible particle 

mass m^. The blue dashed line portrays the red dashed line and the red 

thin dotted lines intersect at the value of true masses. This plot clearly illustrates that because of 
the on-shell constraint, M 2 Cons attains a larger value, even bigger than the corresponding 
once is different from the true mass m^. The most compelling observation about this plot is the 
appearance of a kink exactly at the true mass point for which can be used solely for 

measuring both mp and simultaneously. The lower plot describes the difference - 

with respect to the trial invisible particle mass m^. As expected, the difference between 
both the end points is zero at the true invisible particle mass. 

Although there is no analytic formula in support of the above empirical observation, we 
verified it by checking the slope numerically before and after the true mass point. The 
presence of this kink is also authenticated by various mass spectrums. Also note that unlike 
the constrained variable cannot increase forever with the increase 

of the trial invisible particle mass rh^, owing to the additional heavy resonance mass-shell 
constraint. can maximally reach up to half of the resonant mass and after that it 

would be unphysical. We have not studied the effects of ISR or UTM on this kink solution, 

constraints extend/shift their distribution end-point value over and above the end point where this relative 
constraint is absent. Similar to our case, this can happen when trial mass deviates from true the invisible 
particle mass. Hence, one expects formation or consolidation of similar kink structure. 


rrip = 200 GeV 

/ 

= 100 GeV 

/ 


.... ••. 
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but one expects that the presence of those extra transverse momenta will sharpen the kink 
structure, leaving these realistic studies for future work. 

To demonstrate this behavior in a more quantitative sense, we once again consider the 
toy process with the antler topology with the aforementioned mass spectrum. The top plot 
of Fig. [^depicts the dependence of both as well as the constrained 

with respect to the trial invisible particle mass m^. The red thin dotted lines showing true 
mass lines intersect at the true mass point, mp} = {100, 200} in GeV. This plot clearly 
illustrates that because of the on-shell constraint, M 2 Cons attains a larger value, even bigger 
than the corresponding once is different from the true mass m^. However, both 

of these maximum quantities attain the same value precisely at the true mass. The most 
compelling observation about this plot is the appearance of a kink exactly at this point for 
which can be used for measuring both masses mp and simultaneously. The 
bottom plot describes the variation of difference between two maximums i.e. - 

with respect to the trial invisible particle mass As expected, this difference 
between both the end points should ideally be zero at the true invisible particle mass. 


IV. MORE ASPECTS OF KINK MEASUREMENT 

One of the signihcant challenges with most of the mass variables is the detection of the 
distribution end point, which can reach to the theoretical maximum only after using a large 
amount of data with signihcant statistics. The problem comes from the fact that negligible 
amount of events typically contribute towards the distribution end point. M 2 Cons is also not 
an exception, forming a tail in the distribution towards its maximum value.® 

This feature is clarihed in Fig. where the density of events is displayed as a percentage 
of total data contributing to the M 2 distributions. As a function of the trial invisible mass 
m^, the left plot shows the reference density for the M 2 distribution which can reach up 
to a maximum, above which there are no events and it remains white. For the same data 
set, the right plot shows the density of events for the constrained variable M 2 Cons- The 
color coding represents the percentage of events per 2 GeV bin in M 2 . These upper end 
points are equivalent to the maximum curve in Fig. in the last section. Gompared to 

® On the contrary, at the true invisible mass, M 2 Cons produces a sharper end point as demonstrated in 
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FIG. 5: Density of events as a percentage of total data contributing to M 2 distributions as a 
function of the trial invisible mass my.. On the left, the reference figure is shown for the M 2 
variable which does not include heavy resonance constraint. The similar figure on the right is for 
the constrained variable M 2 Cons) where constraint refers to the G mass-shell constraint. The color 
coding represents the percentage of events per 2 GeV bin in M 2 . Since this distribution reaches 
a maximum for every trial value of above which there are no events, this disallowed range 
kept as white. This upper end point in each plot represents the maximum curve shown in Fig. 
last section. The presence of the kink can be clearly seen from the figure, and it is solely because 
of the on-shell heavy resonance constraint. But the reconstruction of the kink can be challenging 
due to the much smaller number of events at the endpoint, specifically when away from the kink. 
Also, it is interesting to note the changes in event density due to the application of an additional 
constraint. Evidently, a significant number of events shifted towards the end point at the true 
mass, as can be observed in the figure. 

the left figure, M 2 Cons developed a clear kink solely because of the on-shell constraint of 
the heavy resonance. One notices that a tiny fraction of events is actually contributing at 
the end point, specihcally when away from the kink position. Also, it is interesting to note 
the changes in event density due to the additional constraint. At the true mass (kink), a 
signihcant number of events shifted towards the end point, as is also observed in Fig. This 
demonstration is also generated considering the toy process with the antler topology with 
the aforementioned mass spectrum {rriG^mp^m^} = {1000,200,100} in GeV. 
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FIG. 6: The effectiveness of the variable n{m^) as defined in Eq. 11 is in identifying the minimum 


and thus measuring the true mass of the invisible daughter. The function as a percentage of the 
event fraction clearly shows a sharp minimum at rhy. = rriy^. So by identifying the minimum of 
n{my^), one can measure the invisible particle mass accurately. The red band shows the error 
accounting only for the statistical uncertainty. 


We pointed out and discussed the difficulty with determining the end points, which is in 
no way a shortcoming for this variable only. Fortunately, in this present case, the ability to 
simultaneously identify both M 2 and M 2 Cons provides a solution for effectively pointing out 
the kink using all the events, not just relying on the events at the maximum. 

We have already discussed in Sec. [IT] that the additional constraint pushes the M 2 Cons 
towards the higher value compared to M 2 , such that, as long as the trial invisible mass 
is unequal to the true mass there can be enough events generating a larger M 2 Cons than 


Moreover, it is clear from Eq. 10 that the coincides with at the true 


invisible mass. This enables us to dehne a dimensionless variable pointing out the position 
of kink, in a way that was originally proposed in |32]. For a given rh^, one counts all the 
events having M 2 Cons value larger than the corresponding to get the fraction, 

1 1 ^ 

n(m^) = — 'Hi{,M 2 Cons{.'m^) - M^'^'"{m^)). (11) 


i=l 


Here i is the event index with total number N. M{m^ is the number of events in which 


17 










M2Cons{^x) > for any given satisfied by the Heaviside step function, 


My) = 


0, a y<0 
1, if y> 0. 


( 12 ) 


It is easy to follow from Eq. 10 that the quantity n{m^) should ideally be zero at the true 
mass since both M 2 Cons and M 2 share the same maximum value at that point. However, 
on both sides away from this point, substantial events contribute above the hence, 

n(m^) poses a sharp minimum at the true mass point. Considering other realistic effects, 
such as backgrounds, mass width, and experimental errors, can lift the minimum from zero. 
These effects are not considered in this present analysis. However, it is safe to assume that 
the position of the functional minimum can be correctly identihed to get the true invisible 
mass. The advantage of using n{m^) is that it does not rely on some isolated event at 
the end point but rather it relies on a signihcant number of events distributed on a band 
in a two-dimensional plane between and which contribute to establish this 

minimum. 

In our example, theoretical prediction of the function n{fh^) (as a fraction of total events) 
is shown in Fig. The red band is the error accounting only for the statistical uncertainty. 
One can clearly identify the minimum and justify the relation 


n{m^ = m^) = n^inijn^) (13) 

to measure the invisible particle mass accurately. Hence, it is straightforward to measure 
both the parent and daughter masses simultaneously. 


V. RECONSTRUCTION CAPABILITY OF EVENTS 

In this section, we are inclined to explore the event reconstruction capability coming 
from the constrained mass variable M 2 Cons, typically once the invisible particle mass is 
determined, as in last section. Event reconstruction is extremely important in the case of 
spin, polarization and the coupling determination of new physics as well as SM processes with 
the Higgs and top. However, it is almost impossible to determine them exactly for a scenario 
involving multiple invisible particles in a hadron collider, especially for a topology with a 
short decay chain. Attempts have been made to reconstruct events using the transverse 
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FIG. 7: Capability of reconstructing missing daughters momenta is demonstrated using con¬ 
strained variable. Left figure is displaying a normalised distribution of , i with Agi^ = 

Lit I 

~ hence parameterize the deviation from true momenta for the transverse part, With 

”t” either x or y-component of momenta. This reconstruction of momenta is done from the min¬ 
imisation of the M 2 Cons with the true mass of the invisible particle as input. The reconstructed 
invisible momenta is unique and very well correlated with true momenta as the distribution of 

. I has a sharp peak at zero. In a process where the invisible particle mass is unknown, n{m^) 

Lit I 

can be used for invisible particle mass determination and then event reconstruction using M 2 Cons- 

Similarly, right figure displays a normalised distribution of , i with Agi^ = for more 

Liz I 

troublesome longitudinal momentum. Once again, qi^ reconstruction is unique and well correlated 
with true longitudinal momenta of invisible particle. 

mass variable Mt 2 [591 - I6T] known as the M'r 2 -assisted on-shell (MAOS) method in which 
the transverse momenta of the invisible particle are determined from the minimization of 
Mt 2 , and the longitndinal components are determined by solving mass-shell the constraints 
of parent and daughter. In Ref. |17l HE], it is shown that Smin and its constrained sisters 
and can also be used for event reconstruction especially in antler topology, where 
constraints refer to missing transverse momenta and available mass-shell constraints in an 
event. The reconstructed momenta of the invisible particles are derived at the extremum of 
s and the constrained s variables. 

In this section we reconstruct all components of the invisible particle momenta from 
minimization of the (1 -|- 3)-dimensional variable M 2 Cons with the true mass of the invisi¬ 
ble particle as input. The capability of reconstructing the missing daughters momenta is 
demonstrated in Fig. [fusing this constrained variable. The left plot displays a normalized 
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distribution of the deviation (or error) in this reconstructed quantity from that of the true 
value, in the case of transverse part of the invisible momenta. This deviation is parametrized 
using a ratio dehned as 


^reco 


^True 

Hit 


\<i 


True I 

It I 


(14) 


where the snbscript ”t” refers to the transverse (x or y) component of momenta. The re- 
constrncted invisible momenta are proved to be nniqne and very well correlated with the 
trne momenta, as the observed distribntion has a sharp peak at its true value (he. zero 
deviation). Similarly, the right hgure displays a normalized distribntion of the correspond¬ 
ing variable for longitndinal momentnm, and once again one gets a nnique reconstrnction 
well correlated with the true longitudinal momenta of the invisible particle. In a process 
where the invisible particle mass is nnknown, n{m^) can be used for invisible particle mass 
determination, and then events can be reconstrncted using M 2 Cons- 


VI. SUMMARY AND CONCLUSIONS 

Looking forward to another breakthrongh in the second phase of its journey, the LHC is 
discovering credible hints for new physics. Many of the popular BSM models are extended 
with massive exotic dark matter particles. As a result, common signatures, coming from 
snch models at a high-energy collider, typically are detectable SM particles along with a 
pair of invisible particles. They are too complex to measnre all the nnknown masses or to 
fnlly reconstruct those events at the large hadron collider. 

In this paper, we stndy one snch class of events prodnced throngh antler topology. They 
are common in SM Higgs and several BSM scenarios, where the heavy resonance is prodnced 
before semi-invisible decay into visible decay prodncts together with invisibles which can 
either be SM nentrinos or some exotic dark matter particles. Onr objective is to determine all 
the nnknown masses, including the dark matter particles produced from the heavy resonance. 
We consider a new constrained variable M 2 Cons extending the (1 -|- 3)-dimensional mass 
variable M 2 , by implementing additional heavy grandparent mass-shell constraint in the 
minimization. 

This new variable M 2 Cons contains several interesting features. We demonstrate how this 
variable acqnires an event wise higher valne owing to this constraint. In particular, we show 
how this variable moves closer to the unknown parent mass. In addition, the calcnlated 
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invisible momenta at this minimum can provide a close estimate of the true momenta of 
the invisible particles for such events. Both these characteristic features are highlighted and 
exploited further to sharpen the measurements. 

Another striking feature comes out once we analyze the distribution maxima of this 
new variable as a function of the trial values of yet unknown dark matter 

particle mass. This is constructed in an analogy with the popular study of 
which gives a useful correlation curve relating the parents mass with the invisible particle 
mass. But now, under mass-shell constraint, develops a new kink solution over 

the correlation curve exactly at the value where this trial mass coincides with the true mass. 
Hence, this opened another new avenue that produces a kink feature to measure both masses 
simultaneously. 

To handle the sparseness of events towards the distribution end point, we analyze with an 
experimentally feasible observable n(m^) by utilizing both constrained and unconstrained 
variables. This observable does not rely on isolated events at the end point, but instead uses 
a signihcant amount of available data to pinpoint the unknown invisible particle mass from 
the sharp minimum. 

Our method provides a complementary procedure to earlier antler studies and is ap¬ 
plicable to any mass region. We demonstrate our analysis in the large mass gap region, 
considered as a difficult region for the kinematic cusp method. In this region, the cusps 
of many variables are not very sharp, which makes the mass determination more prone to 
error. But the present method can be used safely for better accuracy. We also investigate 
the event reconstruction capability of M 2 Cons, and we reconstruct the unknown invisible 
particle momenta at the constrained minimization. The reconstructed momenta are found 
to be unique and well correlated with the true invisible momenta. 
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